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Discovery of the 3.9 days orbital periodicity in Very High Energy 7-rays. 
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Abstract Recent observations of the binary system 
LS 5039 with the High Energy Stereoscopic System 
(H.E.S.S.) revealed that its Very High Energy (VHE) 
7-ray emission is modulated at the 3.9 days orbital pe- 
riod of the system. The bulk of the emission is largely 
confined to half of the orbit, peaking around the inferior 
conjunction epoch of the compact object. The flux modu- 
lation provides the first indication of 7-ray absorption by 
pair production on the intense stellar photon field. This 
implies that the production region size must be not sig- 
nificantly greater than the gamma-gamma photosphere 
size 1 AU), thus excluding the large scale coUimated 
outflows or jets (extending out to ~ 1000 AU). A hard- 
ening of the spectrum is also observed at the same epoch 
between 0.2 and a few TeV which is unexpected under 
a pure absorption scenario and could rather arise from 
variation with phase in the maximum electron energy 
and/or the dominant VHE 7-ray production mechanism. 
This first-time observation of modulated 7-ray emission 
allows precise tests of the acceleration and emission mod- 
els in binary systems. 

Keywords gamma rays: observations • X-rays: bina- 
ries • individual objects: LS 5039 (HESS J126-148) 
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1 Introduction 

In the commonly accepted paradigm, microquasars con- 
sist of a stellar mass black hole fed by a massive star. 
They can exhibit superluminous radio jets[22|, and hints 
for the presence of an accretion disk. These scaled down 
versions of Active Galactic Nuclei (AGN) also show, due 
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to their much lower mass, variability on timescales shorter 
by several order of magnitudes (down to minutes or even 
seconds) that can be used to constrain the accretion and 
ejection scenarios (e.g. ^^). 

It has been suspected for a long time that these ob- 
ject could emit, through similar acceleration mechanisms 
as in AGN (leptonic or hadronic processes) high energy 
radiation up to the Very High Energy gamma-ray (VHE) 
domain, and that VHE radiation could give insight on 
the very central engine. The discovery of VHE emis- 
sion from LS 5039 by HESS|1] and shortly after from 
LS I +61° 303 by MAGICS' confirmed this longstand- 
ing issue and established VHE astronomy as a powerful 
diagnostic probe of these objects. 

LS 5039, identified in 1997[5B| as a massive X-ray bi- 
nary system with faint radio emission |19J. was resolved 
in 2000 27 into bipolar mildly relativistic radio jets {v ^ 
0.2 c) emanating from a central core, thus placing it 
into the microquasar class. The detection of radio and 
variable X-ray emission TO' and its possible association 
with the EGRET source 3EG J1824-1514,2S suggests 
the presence of multi-GeV particles possibly accelerated 
in jets. LS 5039 is indeed the only object simultaneously 
detected in X-ray and radio in the field of view of the 
unidentified EGRET source |5T]. 

The binary system LS 5039 (Fig^l consists of a mas- 
sive 06.5V star in a ~ 3.9 day mildly eccentric orbit 
(e = 0.35) II around a compact object whose exact 
nature (black hole or neutron star) is still under de- 
bate. Under the assumption of pseudo-synchronization, 
Casares and collaborators 11 constrain the compact ob- 
ject mass in the black hole range {Mx — ^.T^ii'^MQ) and 
obtain a low system inclination {i ^ 25°), but a neutron 
star at higher inclination (i ~ 60°) might still remain 
possible (14,. 



2 H.E.S.S. Observations 

The High Energy Stereoscopic System (H.E.S.S.) is an 
array of four identical Atmospheric Cherenkov Telescopes 
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Fig. 1 Orbital geometry of tiie binary system LS 5039 
vi ewe d from above and using tlie orbital parameters derived 
in|ll|. From Sliown are: pliases (0) of minimum (perias- 
tron) and maximum (apastron) binary separation; epocli of 
superior and inferior conjunctions occurring when tlie com- 
pact object and tiie star are aligned along the observer light- 
of-sight. 

(ACT) 01 and references therein] located in the Southern 
Hemisphere (Namibia, 1800 m a.s.l.) and sensitive to 7- 
rays above 100 GeV. 

LS 5039 was serendipitously detected in 2004 during 
the H.E.S.S. galactic scan(2j. The 2004 observations have 
been followed up by a deeper observation campaign[3] in 
2005, leading to a total dataset of 69.2 hours of observa- 
tion after quality selection. To optimize the coverage over 
the orbit, the observations were spread over more than 
six months, resulting in a wide range of observation con- 
ditions. The observation zenith angles are in particular 
distributed between 5° and ^ 65°, resulting in a trig- 
ger threshold varying between 100 GeV and ~ 1 TeV. 

Data were analysed using two separate calibrations ^ 
and analysis pipelines. The results presented here are 
based on the log-likelihood comparison of the shower im- 
ages with a precalculated semi-analytical mo del [T^. 



3 Results 

After selection cuts, a total of 1969 7-ray events were 
found within 0.1° of the VLBA radio position 
of LS 5039 30 , leading to a statistical significance of 40cr 
(Fig. 13). The best fit position is, in Galactic Coordinates, 
/ = 16.879°, b — —1.285° with statistical and systematic 
uncertainties of respectively ±12" and ±20". This po- 
sition is compatible within la with the VLBA position 
(denoted as a blue star in Fig. ^ and with the Chandra 
source. We obtain an upper limit of 28" (at la) on the 
VHE source extension. 



■D 



-0.5 



-1 



-1.5 



-2 

18 17.5 17 16.5 16 

I (deg) 

Fig. 2 H.E.S.S. excess sky map around LS 5039, smoothed 
by the instrument point spread function. The blue star de- 
notes the position of the VLBA source. The yellow contours 
correspond to the 68%, 95% and 99% confidence level region 
of the EGRET source 3EG J1824-1514. The extended source 
HESS J1825-137 observed in the same field of view can serve 
as a cross-check for timing-analysis. 

3.1 Timing Analysis 

The runwise VHE 7-ray flux at energies > 1 TeV was de- 
composed into its frequency components using the Lomb- 
Scargle periodogram|3S] (Fig. 1) which is appropriate 
for unevenly sampled datasets such as those collected 
by H.E.S.S. 

In order to reduce the effect of the varying instrument 
threshold, all events were used in the lightcurve determi- 
nation and the runwise flux normalisation was extracted 
under the assumption of an average photon index derived 
from all data (F = 2.23 for dN/dU cx E'^. The average 
index assumption in this method notably increases the 
statistics and contributes only to a small error on the 
derived flux above 1 TcV). 

An obvious peak in the Lomb-Scargle periodogram 
occurs at the period 3.9078 ± 0.0015 days, consistent 
with the orbital period determined by Caseres et aL|ll| 
(3.90603 ± 0.00017) and excluding the earher period of 
4.4267±0.0005 days determined by McSwain et all20r21!. 
The peak is highly significant, with a chance probabil- 
ity of ^ 10"^" before trials (estimated via Monte-Carlo 
simulation of random fluxes time-series at the observa- 
tions times as well as random shuffling of observed fluxes) 
and less than 10~^^ after trials. Other peaks with chance 
probability less than 10^''-10~* are present in the peri- 
odogram. 

FigO middle panel, shows the effect of subtracting 
a pure sinusoid at the orbital period. The orbital peak 
disappears as expected, but also the numerous satellite 
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Fig. 3 Lomb-Scargle (LS) periodogram of the VHE runwise 
flux of LS 5039 above 1 TeV (Chance probabihty to obtain 
the LS power vs. frequency). From Zoom: inset around 
the highest peak, which corresponds to a period of 3.9078 ± 
0.0015 days. This period is found to be co mpat ible with the 
orbital period determined by Caseres et al.\ll\ and denoted 
as a red line on the inset. Middle: LS periodogram of the 
same data after subtraction of a pure sinusoidal component 
at the orbital period of 3.90603 days. The orbital peak is 
removed as expected, but also the satellite due to beat of the 
orbital period with various instrumental periods (see text). 
Bottom: LS periodogram obtained on HESS J1825-137 using 
the same dataset (HESS J1825-137 is observed in the same 
field of view.) 
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Fig. 4 Phasogram (Integral run-by-run 7-ray flux above 1 
TeV as function of orbital phase) of LS 5039 from H.E.S.S. 
data from 2004 to 2005, using the orbital ephemeris of llll . 
Each run is ~ 28 minutes. Two full phase periods are shown 
for clarity. The vertical blue arrows depict the respective 
phases of minimum {periastron) and maximum (apastron) bi- 
nary separation. The vertical dashed red lines show the re- 
spective phases of inferior and superior conjunction, when the 
star and the compact object are aligned along the observer's 
line of sight. From 



3.3 Spectral Modulation 



peaks, thus confirming that these peaks are beat periods 
of the orbital period with the various gaps present in the 
H.E.S.S. dataset (day-night cycle, moon period, annual 
period). The bottom panel of the same figure shows the 
result obtained on the neighbouring source HESS J1825- 
137 observed in the same field of view as LS 5039. The 
HESS J1825-137 periodogram doesn't show any statisti- 
cally significant peak, thus demonstrating that the ob- 
served periodicity is genuinely associated with LS 5039. 



3.2 Flux Modulation 

The runwise Phasogram (Fig^j) of integral flux at ener- 
gies > 1 TeV vs. orbital phase {(j>) shows an almost sinu- 
soidal behaviour, with the bulk of the emission largely 
confined in a phase interval (j) ~ 0.45 to 0.9, covering 
about half of the orbital period. The thick red line in Fig 
^represents the component at the orbital frequency de- 
termined with the Lomb-Scargle coefficients. The emis- 
sion maximum (0 ~ 0.7) appear to lag behind the apas- 
tron epoch and to align better with the inferior conjunc- 
tion {(j) = 0.716), when the compact object lies in front 
of the massive star (see Fig. ^ . The VHE flux minimum 
occurs at phase {4> ~ 0.2), slightly further along the orbit 
than superior conjunction (0 = 0.058). 

Neither evidence for long-term secular variations in 
the VHE flux independent of the orbital modulation nor 
any other modulation period are found in the presented 
H.E.S.S. data. 



Due to changing environment with orbital phase (mag- 
netic field strength, stellar photon field, relative position 
of compact object and star with respect to observer, 
. . . ) , the VHE 7-ray emission spectrum is expected to 
vary along the orbit. In such a binary system, the com- 
pact object environment can be modeled with a rela- 
tively good accuracy, and the spectral modulation can 
therefore serve as a quite important diagnostic tool for 
disentangling the possible acceleration, cooling and ab- 
sorption processes. 

We first define two broad phase interval: INFC cen- 
tered on the inferior conjunction (0.45 < (f) < 0.9) and 
its complementary SUPC centered on the superior con- 
junction, corresponding respectively to high and low flux 
states. FigOshows the VHE spectral energy distribution 
for these two phase intervals. The high state is consistent 
with a hard power law with index 7^ = 1.85 ± O.OGstat i 
O.lgyst with and exponential cutoff at Eq — 8.7±2.0 TeV. 
In contrast, the spectrum for low state is compatible with 
a relatively steep {F = 2.53 ± O.Q6stat ± O.lsyst) pure 
power law extending from 200 GeV to ~ 20 TeV. The 
spectral shapes of these two states are mutually incom- 
patible at the level of ~ 2 x 10~^. Interestingly, the flux 
appears to be almost unmodulated at 200 GeV as well 
as around 20 TeV, whereas the modulation is maximum 
around a few (~ 5) TeV. 

When going to smaller phase interval (0.1 phase bins), 
the statistics at high energy becomes too low to effi- 
ciently distinguish between power-law and more compli- 
cated shapes. Fig El shows the results (photon index and 
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Fig. 5 Very high energy spectral energy distribution of 
LS 5039 for the two broad orbital phase intervals defines in 
the text: INFC 0.45 < < 0.9 (red circles) and its com- 
plementary SUPC (blue triangles). The shades regions rep- 
resent the la confidence bands on the fitted functions. Both 
spectral are mutually incompatible at the level of ~ 2 x 10~^. 
A clear spectral hardening is occurring in the 200 GeV to a 
few TeV range during the INFC phase interval. From 




Photon Index 

Fig. 7 Differential flux at 1 TeV vs. photon index. The red 
line shows the best linear fit (correlation factor r = 0.8). 

was found to vary between 1.7 and 2.2 with time. How- 
ever, the X-ray phasogram exhibited a different picture 
than the VHE one, with a flux maximum at ((> ~ 0.2 
(close to the VHE flux minimum) and a second peak 
around (f> ~ 0.8 better ahgned with the VHE flux. 
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Fig. 6 Top: Fitted pure power-law photon index vs. phase 
interval of width A(j> = 0.1. Bottom: Differential flux at 1 TeV 
for the same phase interval. From 0. 



differential flux at 1 TeV) of a pure power-law flt of the 
high energy spectra in 0.1 orbital phase bins (restricted 
to energies below 5 TeV to avoid systematic effect intro- 
duced by the high state cutoff). 

The flux normalisation and photon index are strongly 
correlated, the flux being higher when the spectrum is 
harder and vice- versa, as shown in Figd The correlation 
factor is found to be r 0.8. The photon index varies 
by more than one unit along the orbit, whereas the flux 
normalization varies by a factor of more than 5. Inter- 
estingly, a similar effect, however in a smaller variation 
range, was found in X-ravspH] where the photon index 



4 Interpretation 

4.1 Gamma- Ray production 

The basic paradigm of VHE 7-ray production requires 
the presence of particles accelerated to multi-TeV en- 
ergies and a target comprising photons (for 7-ray pro- 
duction through Inverse Compton effect) and/or matter 
of sufficient density (for 7-ray production through pion 
decay in hadronic processes, e.g. Several model 

classes are available to explain VHE emission from mi- 
croquasars, differencing one from the other by the na- 
ture of accelerated particles and/or the location of the 
acceleration region. In jet-based models, particle accel- 
eration could take place directly inside and along the 
jet, e.g. and references therein], and also in the jet 
termination shock regions [TB|. Non-jet scenarios are also 
available, e.g. 18 , 14 , where the emission arises from the 
interaction of a pulsar wind with the stellar companion 
equatorial wind. 

In LS 5039, the strong observed modulation provides 
new information about the physical processes in micro- 
quasars, placing in particular strong constraints on the 
location of the acceleration region. 



4.2 Modulated absorption by pair creation 

VHE 7-rays produced close enough to the stellar com- 
panion will unavoidably suffer severe absorption via pair 
production (e+e") on its intense photon field. Due to 
the angular dependence of the pair creation cross-section 
and threshold, the optical depth will strongly depend on 
the alignment between the 7-ray production region, the 
star and the observer, leading to an orbital modulation 
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Fig. 8 Comparison of the observed phasogram (black points) 
with the expected flux modulatio n u nder a pure absorption 
scenario (blue line, adapted from fH)') and under a accretion 
disk-jet scenario (red line adapted from In the accretion 
disk-jet scenario the modulation of the accretion rate of the 
stellar wind along the orbit is responsible for a modulation of 
the jet injection and therefore for a modulation of the particle 
acceleration in the jet. 




Fig. 9 Expected effect of the absorption on the j-ray spec- 
trum for different orbital phases. Adapted from The blue 
line denotes the primary unabsorbed spectrum. The black 
lines show the absorbed spectra for various orbital phases 
range between 0° and 180° (bottom to top). The sharply 
peaked pair production cross-section is expected to produce 
a dip in the observed 7-spectrum just above the pair produc- 
tion production threshold, in the 200 — 500 GeV. 

of the VHE 7-ray flux !iiO"i>5'"S"TTT. For the orbital ge- 
ometry of LS 5039, a flux minimum is expected at the 
phase of inferior conjunction, where the effect of the min- 
imum absorption threshold and the maximum column 
density adds up. The overall expected modulation for an 
emission close to the compact object (Fig. |H1 blue line 
adapted from agrees quite well with the observed 
picture, suggesting that absorption plays an important 
role in the observed modulation. 

However, its has to be noted that a signal is ob- 
served by H.E.S.S. around phase = 0, which is unex- 
pected under a pure absorption scenario. Detailed treat- 
ment of pair cascades, e.g. |7], in which the VHE 7-ray 
energy is reprocessed towards somewhat lower values, 



might help solving this issue. Another key expectation 
from the absorption scenario is that the strongest ab- 
sorption, and hence modulation, should occur in the en- 
ergy range E ~ 0.2 to 2 TeV depending on orbital 
phase. Fig. El adapted from [H|, shows the expected evo- 
lution of the VHE spectrum with phase. The flux dip at 
absorption maximum will produce a spectral softening 
at low energy ('^< SOOGeV) around phase (/) 0, and a 
spectral hardening above. This picture is not consistent 
with the observed spectral modulation (Fig. I^J with no- 
tably an almost unmodulated flux at ^ 0.2 TeV, thus 
suggesting that additional processes must be considered. 



4.3 Modulated particle acceleration 

VHE 7-ray production can be produced by accelerated 
electrons through the inverse-Compton (IC) scattering 
of stellar photons of the companion star, and/or acceler- 
ated hadrons through their interaction with surrounding 
photons and particles. In this scenario, the efficiency of 
VHE 7-ray production will peak around periastron, re- 
flecting the higher target photon density. 

The high temperature of the companion star means 
that IC 7-ray production proceeds primarily in the deep 
Klein-Nishina regime (where the IC cross-section is sharply 
reduced compared to the Thompson regime). Due to 
changing magnetic field strength and target photon den- 
sity along the orbit, the maximum electron energy will 
be increased by a factor of roughly 10 between periastron 
and apastron, where the IC cooling becomes less efficient 
(for a more detailed discussion, see 3 ). This might agree 
quite well with the observed spectral modulation. More- 
over, synchrotron cooling takes over IC cooling above 
some energy e « 6[(i3/G)(d/i?*)]~^ changing with or- 
bital phase. Much stronger synchrotron losses[23| would 
then produce a phase dependent spectral break, with a 
break energy lower by a factor of ~ 10 at periastron. 

Other effects such as angular dependence of IC scattering [TTI 
could also introduce spectral hardening at apastron phase. 



4.4 Modulated particle injection 

In an accretion disk-jet scenario, orbital modulation could 
arise from modulation of the jet injection rate. Paredes 
et al'SS' propose a leptonic jet model for LS 5039, where 
a slow equatorial stellar wind induces a modulation of 
accretion rate (maximum shortly after periastron). An 
additional stream wind formed next to periastron is re- 
quired to reproduce the doubled-peak X-ray lightcurve. 

Their predicted VHE orbital modulation (Fig|H| red 
line) does not agree well with the observed orbital mod- 
ulation. In particular, the enhancement in accretion rate 
shortly after periastron is not observed. However, com- 
plication arise from the detail treatment of absorption 
effect, which is this case depends on the assumption of 
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the position of 7-ray production along the jet. In par- 
ticular, emission in the small scale jets would produce a 
similar absorption pattern as in FiglSJ blue line. 



5 Conclusion 

New observations by HESS have established orbital mod- 
ulation of the VHE 7-ray flux and energy spectrum from 
the XRB LS 5039. The flux vs. orbital phase profile pro- 
vides the first indication for 7-ray absorption within an 
astrophysical source, suggesting that a large part of the 
VHE 7-ray production region lies inside the pair absorp- 
tion photosphere (within ^ 1 AU) around the massive 
stellar companion. However, not all of the observed ef- 
fects can be explained by absorption alone. A detailed 
study is now required to fully explain these new observa- 
tions and understand the complex relationship between 
7-ray absorption and production processes within these 
binary systems. 
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